Unified description of charge and spin excitations of stripes in cuprates. 
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Abstract 

We study stripes in cuprates within the one-band and the three-band Hubbard model. Magnetic and charge excitations 
are described within the time-dependent Gutzwiller approximation. A variety of experiments (charge profile from 
resonant soft X-ray scattering, incommensurability vs. doping, optical excitations, magnetic excitations, etc.) are 
described within the same approach. 
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1. Introduction 

If the mechanism of high-T c superconductivity is 
electronic, to understand the excitation spectrum is 
as important as to understand phonons was impor- 
tant to develop BCS theory. In this regard charge 
and spin inhomogeneous states, often found in 
strongly correlated systems, are interesting because 
they can support new collective modes, "electronic 
phonons" , that would not be present in a weakly 
interacting fluid. For example stripes in cuprates 
have many oscillatory modes that can couple with 
carriers and eventually lead to pairing and super- 
conductivity or anomalous Fermi liquid properties. 

Cuprates, unlike BCS superconductors, are 
unique and therefore their remarkable proper- 
ties may well be rooted on details. For example 
a recent study shows that superconductivity in 
La2_ x Sr x Cu04 appears only when incommensurate 
low energy scattering parallel to the CuO bond 
is present and not, when the scattering is on the 
diagonals (at very low doping), or when the sys- 
tem is in a more conventional Fermi liquid phase 
at high doping [1], thus the excitation spectra must 



be understood in detail. In the last years we have 
developed a method to compute collective excita- 
tions of inhomogeneous strongly correlated systems 
based on a time-dependent Gutzwiller approxima- 
tion (GA), named GA+RPA[2,3,4,5,6]. 

Many theoretical descriptions of cuprates are 
suited for only one experiment. Here we show that 
within the same approach one can reproduce a va- 
riety of experimental features of La2- x Sr x Cu04 in 
the charge and magnetic channel. Results in the 
one-band Hubbard model (1BHM) and the three- 
band Hubbard model (3BHM) are quite similar. We 
compare results from both models. 

For the 3BHM we use the following LDA param- 
eter set: e p — ed = 3.3 eV for the splitting between 
the diagonal energies of a hole in the copper d- and 
oxygen p orbital, t p d = 1.5 eV (t pp — 0.6 eV) for 
the p-d (p-p) hopping integral, Ud — 9.4 eV (U p — 
4.7 eV) for the repulsion between two holes on the 
same Cu (O) orbital and U p d = 0.8 eV for the Cu- 
O repulsion[7]. For the 1BHM results reported here 
we use t = 353.7meV for the nearest neighbor hop- 
ping, U = 8t for the effective on-site repulsion, and 
t' = — 0.2i for the next-nearest neighbor hopping. 
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Fig. 1. Hole charge minus undoped charge in the 3BHM as 
a function of atomic position in the direction perpendicular 
to the stripes for various doping x, symmetry (BC or SC) 
and charge periodicity d = 9 (a); and d = 4 (b) and (c). 

2. Ground state properties 

Static properties are computed within the 
Gutzwillcr approximation. Probes which are sen- 
sitive to the Cu and O atomic character need to 
be addressed with the 3BHM. In Fig. 1 we show 
the evolution of the charge profile with doping in 
the 3BHM for metallic stripes parallel to the CuO 
bond. The charge profile in real space has a width 
of about 4 lattice sites [8]. The width of the stripes 
defines two regimes: for low doping (a), the stripes 
do not overlap and therefore are weakly interacting, 
for high doping (b),(c) stripes overlap. 

Recently Abbamonte and collaborators have 
shown that our charge profile predicted in the 
3BHM[8] is in excellent agreement with the profile 
measured with resonant soft X-ray scattering [9]. 

The overall shape of the charge profile in the 
1BHM is similar to the one in the 3BHM (see 
Ref. [10]). One also finds a non overlapping regime 
[d > 4) , and an overlapping regime (d < 4) . 

For large and negative —t'/t (non shown here) one 
finds that checkerboard states are favored[12]. 

In Fig. 2 we show the energy for metallic stripes 
as a function of doping for the 1BHM. In each curve 
the charge periodicity perpendicular to the stripe 
is fixed and takes integer values (from left to right) 
d = 10 — 3. In this computations stripes are bond- 
centered (BC). The energy for site-centered (SC) 
stripes is slightly lower at small doping but the dif- 
ferences decreases with the system size. At optimum 
doping both structures become degenerate. 

A small difference with respect to the 3BHM is 
that in the latter BC stripes are more stable at low 
doping and become quasidegenerate with SC ones 
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Fig. 2. Energy per site as a function of doping for stripes in 
the 1BHM. The curves are labeled by the charge periodicity 
d. For clarity we subtracted the energy of the AFM solution 
and the line nox with fio = — 1.6t. Different choices of fio 
correspond to different choices of the origin of the energy of 
the single particle states and do not change the physics. The 
inset reports the incommensurability as obtained from the 
present calculation (line) compared with experiment [11]. 

at optimum doping. Thus the role of SC and BC is 
interchanged. We expect the three-band results to 
be more reliable in this respect. 

In the inset we plot the incommensurability vs. 
doping taking into account the range of stability of 
each solution and compared with experiments from 
Refs. [11]. The result is a Devil's staircase. Up to x » 
1 /8 the plateaux are short and correspond to a num- 
ber of added holes per unit length along the stripe 
close to v ~ 0.5. As doping increases one jumps from 
one solution to the other and the density of stripes 
increases with doping. This explains the behavior 
of the incommensurability e = l/(2d) = x/(2v) w 
x as seen in neutron scattering experiments in this 
doping range. For x > 1/8 the right branch of the 
d = 4 solution is more stable than the v w 0.5 and 
d = 3 solution due to the overlap effect and one gets 
a wider plateau explaining the saturation of the in- 
commensurability seen in the experiment. 

Some details are slightly different in the 3BHM. 
The e — 1/8 plateau starts at lower doping and 
ends at higher doping improving the agreement with 
experiment [8]. For doping x > 1/8 holes populate 
the stripes and therefore the baseline in Fig. 1 in- 
creases. 

For doping x > 0.195 (x > 0.225) in the one-band 
(three-band) model we find the d = 3 stripe (e = 
1/6 « 0.17) to become the lowest energy solution. 
At this doping a variety of different solutions be- 
comes close in energy and the present saddle-point 
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Fig. 3. Dispersion relation of the low-energy transverse exci- 
tations. We show the experimental result for La2Cu04 after 
Ref. [16] and the GA+RPA result. 

approximation brakes down. Probably around this 
doping a quantum melting of stripes takes place. 

It is also possible that the d = 4 stripe solution 
phase separates with the overdoped Fermi liquid 
skipping the d = 3 solution. This scenario will also 
produce a large e = 1/8 plateau in the incommen- 
surate scattering. In this case long range Coulomb 
effects have to be taken into account [13]. 

One can see from the plot of the energy in Fig. 2 
that the electron chemical potential [i = —dE/dx 
is approximately constant in the low doping regime 
(weakly interacting stripes) whereas it decreases in 
the overlapping stripe regime. The same behavior 
was found in a dynamical mean-field study of the 
one-band Hubbard model [14]. This is in qualitative 
agreement with experiment[15] although we find a 
larger rate of change of [i with doping than in the ex- 
periment. It is also possible that this is due to phase 
separation between d = 4 stripes and the overdoped 
Fermi liquid plus long-range Coulomb effects [13]. 

3. Magnetic excitations 

The dispersion relation for magnetic excitations 
in the insulator computed within GA+RPA applied 
to the 1BHM are shown in Fig. 3 and is in excellent 
agreement with the experiment. The dispersion is 
weakly sensitive to t' jt and it is strongly sensitive 
to U/t. A larger U/t produces a flatter dispersion 
in the (it, 0) - (ir/2,ir/2) direction which does not 
agree with the experiment [4]. 

The collective excitations in the stripe phase 
with the same parameter set have been reported 
in Ref. [5]. One obtains a resonance energy at 65 
meV to be compared with the experimental value 
55 mcV[17]. A slight decrease of U improves the 
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Fig. 4. Optical conductivity labeled by doping, system size 
and, in the case of stripes, interstripe distance. The units 
of conductivity are given by tro = 3.6 X 10 2 (f2cm) — 1 with a 
background dielectric constant ej, = 2.6 The curve labeled 
X = 0.028 corresponds to the single-hole solution. For larger 
dopings the figure is an average over the electric field direc- 
tions parallel and perpendicular to the BC stripes. The inset 
shows the MIR band excluding the Drude component. We 
used a Lorentzian broadening of 0.2eV. 

agreement of the magnetic excitations with experi- 
ment as reported in Ref. [6] but produces a charge 
transfer gap in the insulator which is 10% too low. 
The present parameter set is a compromise between 
the charge and magnetic channel. Details of the 
magnetic response and the doping dependence are 
discussed in an accompanying paper [18]. 

4. Optical conductivity 

The optical conductivity is qualitatively similar 
in the three-band and 1BHM. In Fig. 4 we show the 
result in the 3BHM. The charge transfer (CT) gap 
in the insulator is close to 2eV in good agreement 
with the experiment. In the 1BHM the CT gap is 
also at 2eV. This is mainly determined by our value 
of U = 2.8eV confirming our parameter choice. 

Doping induces a doping dependent mid-infrared 
(MIR) band and Drude weight at zero energy since 
our stripes are metallic. The MIR band is due to 
collective lateral fluctuations of the stripe. These 
lateral fluctuations form a band of excitations as for 
a violin string and the MIR band corresponds to the 
zero momentum component. 

The "string" fluctuation band is massive due to 
pinning of the commensurate stripes to the under- 
laying lattice. However, both the experiment [19] and 
our computations [3] show that the mass decreases 
with doping. This is rooted in the quasi-dcgcncracy 
found among BC and SC stripes and indicates that 
stripes form a floating phase at optimum doping. 
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Table 1 

Selected excitations in meV in the one-band Hubbard model 
within the GA+RPA and in Experiment. Values marked 
with * are extrapolated. 





Momentum 


X 


Theory 


Exp. Ref. 


CT gap 








~ 2000 


~ 2000 [19] 


MIR band 





0.08 


~ 500 


~ 250* [19] 


Magnon 


(7r/2 j7 r/2) 





293 


286 ± 5 [16] 


Magnon 


(7T,0) 





326 


333 ±7* [16] 


Resonance 





1/8 


65 


55 [17] 



The MIR band produces an absorption at low en- 
ergy which explains the failure of the Drude model 
at optimum doping. This indicates that anomalous 
Fermi liquid properties are rooted in the low energy 
excitations of the floating phase. 

In the IBHM the MIR band is at energy ~ 0.5 
eV for doping 0.08 which is higher than in the ex- 
periment (Table 1). Experimentally the band is lo- 
cated at ~ 0.5 eV at very low doping and softens 
with doping faster than what we found. This may be 
due to a finite size effect since our optical conduc- 
tivity computations are done in system sizes much 
smaller than the computations of magnetic proper- 
ties (16 x 4 and 40 x 40 respectively.) 

Table 1 summarizes the energy of selected excita- 
tions in the IBHM with the present parameter set. 
It is remarkable that such a simple model with only 
three parameters is able to provide a reasonable de- 
scription of a variety of excitations in different chan- 
nels plus several ground state properties. 

5. Discussion and Conclusions 

In cuprates true magnetic long-range order is 
sometimes detected in experiment [11] however, 
more often, Bragg peaks are not found. Two possi- 
bilities arise: the system may be close to a quantum 
critical point in the disordered phase (dynamic 
stripes) or the system may be in the ordered side 
of the transition but long-range order is not ob- 
served because of disorder (glassy stripes) . Dynamic 
stripes have quasi long-range order in time and 
space whereas glassy stripes at low temperatures 
have practically long-range order in time but short 
range order in space as in a structural glass. This 
latter situation is clear for example in the magnetic 
channel at low doping, in muon spin-relaxation ex- 
periments where below a certain temperature the 
system develops a static (on the muon time scale) 
magnetic field which is not necessarily accompa- 



nied by Bragg peaks [20], thus the spin rotational 
symmetry is broken without long range order. In 
these quantum or classically disordered cases our 
results apply as long as the energy scale is not too 
low (typically tenths of meV). How this reflects on 
the Fermi surface is discussed separately [21]. 

To conclude we have shown that the GA+RPA 
approximation allows for a unified description of col- 
lective modes in the charge and magnetic channel of 
striped cuprates. Results are similar in the one-band 
and the 3BHM although some details may differ in 
which case the 3BHM is expected to be more accu- 
rate. These modes are likely to play an important 
role in the anomalous properties of cuprates. 
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